On polynomial solutions of Heun equation 
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By making use of a recently developed method to solve linear differential equa- 
tions of arbitrary order, we find a wide class of polynomial solutions to the Heun 
equation. We construct the series solution to the Heun equation before identifying 
the polynomial solutions. The Heun equation extended by the addition of a term, 
—a/x, is also amenable for polynomial solutions. 

Heun equation, a second order linear differential equation having four regular singularities, 
is well-known in mathematical literature [1, 2] and has appeared in a number of physical 
problems, like quasi-exactly solvable systems [3], higher dimensional correlated systems [4], 
Kerr-de Sitter black holes [5], Calogero-Moser-Sutherland systems [6], finite lattice Bethe 
ansatz systems [7] etc. We analyze this equation by using a newly developed method to 
solve linear differential equations [8]. This approach to Heun equation is very transparent 
in obtaining solutions; in particular the polynomial solutions come out very naturally as is 
shown below. 

A single variable differential equation, after appropriate manipulation, can be cast in the 
form 

(F(D) + P(x, d/dx))y(x) — , (1) 

where, D = x-^, and F(D) = J2 n a nF> n is a diagonal operator in the space of monomials, 
with a n 's being some parameters. P(x, d/dx) can be an arbitrary polynomial function of x 
and d/dx, excluding the diagonal operators. It can be shown that the following ansatz, 



y(x) 



E(-i) 



m=0 



1 



F(D) 



P(x, d/dx) 



x 



(2) 



is a solution of the above equation, provided, F(D)x x = and the coefficient of x x in y(x)—x x 
is zero, in order to ensure that the solution, y(x), is non-singular [8]. The fact that D is 
diagonal in the space of monomials, x n , makes 1/F(D) well defined in the above expression. 
For explicating the working of the present method, consider the Jacobi differential equation 

[9] 



a - (a + (3 + 2)x]-f- + n{n + a + (3 + 1) 
dx 



p( a >V(x) = 



(3) 



where, P^ Q,/3 ^(x)'s are the Jacobi polynomials. Rewriting the above equation in the form 
of Eq. (1), one finds F(D) = D 2 + (a + (3 + l)D - n(n + a + (3 + 1) and P(x, d/dx) = 
~lx s ~ ~ a )Jr' F(D)x x = yields two solutions A = n and A = —(n + a + (3 + 1). From 
Eq. (2), it is clear that, A = n results in a polynomial solution, 
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(D-n)(D + a + (3 + l) 
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(4) 
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Notice that the choice of A = — (n + a + (3 + 1) results in the other linearly independent 
solution with negative powers of x. We now proceed to study the Heun equation by using 
the above technique. The Heun equation, with four sigularities at x — 0, 1, c and oo, is given 

by 

d 2 ( i 5 e \ d a[3x — q 



+ 



+ 



y{x) = 



dx 2 \x x — 1 x — cjdx x(x — l)(x — c) 
Rewriting the above equation after multiplying it by x, one has 



d 2 d d 2 d 
cx -—7 + 7cx- ( 1 + c)x —t; - (1 + cn + dc + e\x qx 

dx 1 dx dx 1 dx 



d 2 d 
+x 4 — + (7 + 5 + e)a: 3 — + apx 2 



y(x) = 



(5) 



(6) 



In the above equation, 



d 2 

r.r— + -fcx^- = F(D) = cD{D + 7 - 1) 



dx 2 



d 

dx 



and 
where, 



P(x, d/dx) = A+i + A + 2 , 



A +1 = -(1 + c)x z ^r - [(1 + 0)7 + 5c + e}x 2 4~ ~ qx 
dx z dx 



is a degree +1 operator, and 

d 2 d 

A +2 = +^ 4 -^ + (7 + 5 + e ) x3 ^: + a f 3x2 ' 

is a degree +2 operator, respectively, d, the degree of an operator, O, is defined from 
[D, O] = dO. Now, 

F(D)x x = , 

yields A = or 1 — 7. The solutions, for the values of A = and A = 1 — 7, can be written 
as _ 

-(A +1 + A +2 ) 



m=0 



D(D + 7- 1) 



x 



(7) 



The above is an infinite series solution to the Heun equation. Below, we look for the 
polynomial solutions. 

Case (i): Rewriting the Heun equation after dividing it by x, we get, 



D 2 + (7 + 5 + e - 1)D + a(3 - (1 + c)x-^ - [(1 + c) 7 + 5c + e]^- - | 



d 2 Id 

dx 2 ^ x dx 



y{x) = . 



(8) 



Here, 



F(D) = D 2 + (-f + 5 + e-l)D + ap , 
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and 
where, 

and 



P = A_! + A. 



2 , 



d 2 da 
A-! = -(l + c)x— - [(l + c) 7 + 5c + e]— - i , 

ax 2 ax x 

2 ° dx 2 ^°xdx 

The condition, F(D)x x = 0, after imposing the constraint 7 + <5 + e = a + (3 + 1, gives 
A = —a or — /3. In order to have polynomial solutions, one needs to impose q = in 
A_i and furthermore, either —a or — (3 must be a positive integer, say n. For the sake of 
illustration, we choose —a = n. Then, the polynomial solutions to the Heun equation, y n (x), 
can be written as, 



y n (x) = 



k=0 



(D-n)(D + f3) 



(A_ 1 + A_ 2 ) 



x 



(9) 



Case(ii): By performing a similarity transformation on Eq.(7) by x^ 1 ^ and writing y[x) 



x 



1"7, 



> n (x), the choice of q = (5c + e)(7 — 1) yields, 



D 2 + (1 - 7 + 5 + e)D + (1 - 7 ) (5 + e) + a(3 
-(l + c)x^-[(l + c)(2-7)+5c + e]^ 



d" 2 (2 - 7 )c d ' 
+c h 

dx 2 x dx 



here, 
and 



F(D) = D 2 + (1 - 7 + 5 + e)D + (1 - 7 )(<J + e) + a/5 



(10) 



P(x, d dx) = - 1 + c x— - 1 +c 2-7 +& + £- + c— + ^ ^ — 

ax 2 ax dx 1 x dx 



F(D)x x = results inA = 7 — 1 — /3 or 7 — 1 — a, after imposing the constraint 7 + <5 + e = 
a + /3 + 1. In order to have polynomial solutions, either 7 — 1— (3 or 7 — 1 — a must be a 
positive integer. Let 7 — 1 — (3 = n, a positive integer, then the solution is 



m=0 



(D-n)(D-7 + l + a) 



P(x, d/dx) 



x" 



(11) 



It is interesting to note that if one extends the Heun equation by adding a term, —cr/x, 
the above analysis holds and yields polynomial solutions. Starting from 



d 2 f j 5 e \ d afix — q — a/x 

dx 2 \x x — 1 x — cjdx x(x — l)(x — c) 



Y n (x) = , 



(12) 
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and performing a similarity transformation by writing Y n (x) = x 1 ~ J+a Xn{%), and by impos- 
ing the constraint q = (1— r )+a)[{l+c)a+5c+e] (instead of the usual one 7+5+e = a+/3+l), 
one finds that Xn( x ) obeys the following differential equation, 

D 2 + (2a + 1 - 7 + 5 + e)D + (1 - 7 + a)(a + 5 + e) + a(3 

-(1 + c)x^ - {[2(1 + a) - 7 ](1 + c) + 5c + e}^- 
d 2 [2(1 + a) - 7]c d 1 , , 

In the above case, 

F(D) = D 2 + (2a + l- 1 + 5 + e)D + (l- 1 + a)(a + 5 + e) + a(3 , 

and the rest of the terms belong to P(x, d/dx). Imposing F(D)x x = yields, 

A 2 + (2cr + 1 - 7 + 5 + e)A + (1 - 7 + cr)(<7 + 5 + e) + a(3 = , 

which has two roots, A±. In order to have polynomial solutions, either A + or A_ or both 
must be a positive integer. The solutions can then be written as 



x 



E(- 1 )' 



m=0 



(D-\ + )(D-\_) 



P(x, d/dx) 



x 



(14) 



In conclusion, we analyzed the Heun equation by making use of a recently developed 
method to solve linear differential equations of arbitrary order. Apart from the series solu- 
tion, we have obtained polynomial solutions to the Heun equation. Further, we extended 
the Heun equation by adding a term, —a/x, and showed that polynomial solutions also exist 
in this case. The solutions found here may throw new light on physical problems involving 
Heun equation. 
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